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Abstract — Stochastic network calculus is a theory for 
stochastic service guarantee analysis of computer com- 
munication networks. In the current stochastic network 
calculus literature, its traffic and server models are typ- 
ically defined based on the cumulative amount of traffic 
and cumulative amount of service respectively. However, 
there are network scenarios where the applicability of such 
models is limited, and hence new ways of modeling traffic 
and service are needed to address this limitation. This pa- 
per presents time-domain models and results for stochastic 
network calculus. Particularly, we define traffic models, 
which are defined based on probabilistic lower-bounds on 
cumulative packet inter-arrival time, and server models, 
which are defined based on probabilistic upper-bounds 
on cumulative packet service time. In addition, examples 
demonstrating the use of the proposed time-domain models 
are provided. On the basis of the proposed models, the five 
basic properties of stochastic network calculus are also 
proved, which implies broad applicability of the proposed 
time-domain approach. 

I. Introduction 

Stochastic network calculus is a theory dealing with 
queueing systems found in computer communication 
■networks 

[3] [4]. It is particularly useful for analyzing net- 
works where service guarantees are provided stochas- 
tically. Such networks include wireless networks, multi- 
access networks and multimedia networks where appli- 
cations can tolerate some certain violation of the desired 
performance 0. 

Stochastic network calculus is based on properly de- 
fined traffic models [6][3][4][7][8][9] and server models 
00. In the existing models of stochastic network 
calculus, an arrival process and a service process are typ- 
ically modeled by some stochastic arrival curve, which 
probabilistically upper-bounds the cumulative amount 
of arrival, and respectively by some stochastic ser- 
vice curve, which probabilistically lower-bounds the 



cumulative amount of service. In this paper, we call 
such models space-domain models. Based on the space- 
domain traffic and server models, a lot of results have 
been derived for stochastic network calculus. Among the 
others, the most fundamental ones are the five basic 
properties O: (P.l) Service Guarantees including 
delay bound and backlog bound; (P.2) Output Charac- 
terization; (P.3) Concatenation Property; (P.4) Leftover 
Service; (P5) Superposition Property. Examples demon- 
strating the necessity of having these basic properties 
and their use can be found 0. 

Nevertheless, there are still many open research chal- 
lenges for stochastic network calculus, and a critical 
one is time-domain modeling and analysis 0. Time- 
domain modeling for service guarantee analysis has 
its root from the deterministic Guaranteed Rate (GR) 
server model flOl . where service guarantee is captured 
by comparing with a (deterministic) virtual time func- 
tion in the time-domain. This time-domain model has 
been extended to design aggregate-scheduling networks 
to support per-flow (deterministic) service guarantees 
ifTTTl iTTSTl . while few such results are available from space- 
domain models. Other network scenarios where time- 
domain modeling may be preferable include wireless 
networks and multi-access networks. 

In wireless networks, the varying link condition may 
cause failed transmission when the link is in 'bad' con- 
dition. The sender may hold until the link condition 
becomes 'good' or re-transmit. For such cases, it is 
difficult to directly find the stochastic service curve in 
the space-domain because we need to characterize the 
stochastic nature of the impaired service caused by the 
'bad' link condition. A possible way is that we use 
an impairment process to characterize the impaired 
service. However, how to define and find the impairment 
process arises another difficulty. Even though we can 
define an impairment process, we may first convert 



the impairment process into some existing stochastic 
network calculus models, and then further analyze the 
performance bounds. The obtained performance bounds 
may become loose because of such conversion. If we 
characterize the serivce process in the time-domain, we 
can use random variables to represent the time intervals 
when the link is in 'bad' condition. Analyzing the 
stochastic nature of such random variables would be 
easier. In addition, this way can avoid the difference 
introduced by the intermediate conversion. 

In contention-based multi-access networks, backoff 
schemes are often employed to reduce collision occuring. 
Because the backoff process is characterized by backoff 
windows which may vary with the different backoff 
stages, it is quite cumbersome for a space-domain server 
model to characterize the service process with the con- 
sideration of the backoff process. This also prompts the 
possibility of characterzing the service process in the 
time-domain. Having said this, however, how to define a 
stochastic version of the virtual time function and how 
to perform the corresponding analysis are yet open O. 

The objective of this paper is to define traffic models 
and server models in the time-domain and derive the 
corresponding five basic properties for stochastic net- 
work calculus. Particularly, we define traffic models that 
are based on probabilistic lower bounds on cumulative 
packet inter-arrival time. Also, we define server models 
that are based on some virtual time function and proba- 
bilistic upper bounds on cumulative packet service time. 
In addition, we establish relationships among the pro- 
posed time-domain models, and the mappings between 
the proposed time-domain models and the existing space- 
domain models. Furthermore, we prove the five basic 
properties based on the proposed time-domain models. 

The remainder is structured as follows. Sec. [TT] in- 
troduces the mathematical background and fundamental 
space-domain models and relevant results of stochastic 
network calculus. In Sec. Jill we first introduce the 
time-domain deterministic traffic and server models, and 
then extend them to stochastic versions. In addition, the 
relationships among them as well as with some existing 
space-domain models are established. Sec. [TV] explores 
the five basic properties. Sec. Ivl summarizes the work. 

II. Notation and Background 

To ease expression, we assume networks with fixed 
unit lengthQ packets. By convention, we assume that 

'The results can also be extended to networks with variable-length 
packets while the expression and results will be more complicated. 



a packet is considered to be received by a network 
element when and only when its last bit has arrived to 
the network element, and a packet is considered out of 
a network element when and only when its last bit has 
been transmitted by the network element. A packet can 
be served only when its last bit has arrived. All queues 
are assumed to be empty at time 0. Packets within a flow 
are served in the first-in-first-out (FIFO) order. 

A. Notation 

Let p n , r(n), a(n) and d(n) (n = 0,1,2,...) denote 
the n th packet of a flow, its allocated service rate, its 
arrival time and its departure time, respectively. Let A(t) 
and A*(t) respectively denote the number of cumulative 
arrival packets and the number of cumulative departure 
packets by time t. By convention, we assume a(0) = 0, 
d(0) = 0, .4(0) = and A*(0) = 0. For any < s < t, 
we denote A(s, t) = A(t)-A{s) and4*(s,t) = A*(t)- 
A*(s). 

In this paper, a(n) and A(t) will be used to represent 
an arrival process interchangeably. A departure process 
will be represented by d(n) and A*(t) interchangeably. 

The following function sets are often used in this 
paper. Specifically, we use Q to denote the set of non- 
negative wide-sense increasing functions as follows: 

9 = {<?(•) : V0 < x < y, < g(x) < g(y)} 

We denote by Q the set of non-negative wide-sense 
decreasing functions: 

Q = {<?(•) : V0 < x < y, < g{y) < g(x)} 

Let T denote the set of functions in Q, where for each 
function /(•) € T, its nth-fold integration, denoted by 

f n \x) = (J™dy) n f{y), is bounded for \fx > and 
still belongs to T for Vre > 0, or 

T ={/(•): Vn > 0, ( / dy) n f(y)eF}. 
For ease of exposition, we adopt 

[x] + = max[0, x] and [x]\ = min[l,x], 

and assume that for any bounding function f(x), f(x) = 
1 for Vx < 0. 

B. Max-plus and Min-plus Algebra Basics 

An essential idea of (stochastic) network calculus is to 
use alternate algebras particularly the min-plus algebra 
and max-plus algebra lfT3ll to transform complex non- 
linear network systems into analytically tractable linear 
systems JH. To the best of our knowledge, the existing 



models and results of stochastic network calculus are 
mainly under the space-domain and based on min-plus 
algebra that has basic operations particularly suitable 
for characterizing cumulative arrival and cumulative ser- 
vice. For characterizing arrival and service processes in 
the time-domain, interestingly, the max-plus algebra has 
basic operations that well suit the need. 

In this paper, the following max-plus and min-plus 
operations are often used: 

• Max-Plus Convolution of g\ and g 2 is 

{gi®92){x) = sup {gi(y) + g2(x - y)} 

0<y<x 

• Max-Plus Deconvolution of g\ and g 2 is 

(9i&92)(x) = inf {gi(x + y)- g 2 (y)} 
y>o 

• Min-Plus Convolution of g\ and g 2 is 

(gi ® 92)(x) = inf {gi(y) + g 2 (x-y)} 

0<y<x 

• Min-Plus Deconvolution of g\ and g 2 is 

(gi 92)(x) = sup{5i(x + y) - g 2 (y)} 
y>o 

In this paper, when applying supremum and infimum, 
they may be interpreted as maximum and minimum 
whenever appropriate, respectively. 

C. Preliminaries 

The following lemma is often used for later analysis 
and thus listed: 

Lemma 1. For the sum of a collection of random 
variables Z = J27=i no mat ter whether they are 
independent or not, there holds for the complementary 
cumulative distribution function (CCDF) of Z: (See 
Lemma 1.5 in /j^j/) 



F z (z) < F Xl ® ■ ■ ■ ® F Xn (z) 
where Fz = P{Z > z}, —00 < z < 00. 



(1) 



For later analysis, we need some transformation be- 
tween the number of cumulative arrival packets by time 
t, i.e., A(t), and the time of a packet arriving to the 
system, i.e., a(n). 

If A(t) is upper-bounded with respect to some func- 
tion a(t) £ Q, we have the following lemma. 

Lemma 2. For function a(t) £ Q, there holds: 
1) the following statements are equivalent: 

a) VO < s < t, A(s, t) < a(t - s) + x for 
Vx > 0; 



b) V* > 0, A(t) <A® a{t) + x for Vx > 0; 
2) if Vi, x > 0, A(t) < A a(t) + x holds, then we 
have a{n) > a®X{n) — y, where A(n) £ Q is the 
inverse function of a(t) and defined as follows 



A(n) = inf{r : a{r) > n} 



and 



y = sup[A(/c) — A(/c — x)\. 

k>0 



(2) 



(3) 



Proof: (1) For (a) — > (b), from the condition, we 
obtain A(s, t) - a(t - s) - x < for V < s < t. Then, 
there holds 

sup [A(s, t) — a(t — s) — x] < 

0<s<t 

which implies 

A(t) - inf L4(s) + a(t -s)]-x<0. 

0<s<t 

Thus, we conclude A(t) < A <8> a(t) + x for Vt, x > 0. 
For (b) — > (a), from the condition, we have 

A(t) - inf \A(s) + a(t - s)] - x < 

0<s<t 

which implies 

sup [A(s, t) — a(t — s) — x] < 0. 

0<s<t 

Then there must be A(s, t) — a(t — s) — x < for 
V0 < s < t. Thus, A(s, t) < a(t - s) + x holds for 
V0 < s < t and Vx > 0. 

(2) From (1), we know that A(t) < A (8) a(t) + x is 
equivalent to A(s, t) < a(t — s) +x for V0 < s < t and 
Vx > 0. Then for V0 < m < n, we have 

A(a(m), a + (n)) < a{a + {n) — a(m)) + x 

where a + {n) = a(n) + e with e — > 0. We also know 

n — m < A(a(m), a + (n)) < a(a+(n) — a(m)) + x 

Taking the inverse function of a(a + (n) — o(m)) yields 

a+(n)— a(m) > A(n— m— x) 

= A(n — m) — [A(n — m) — A(n — m — x)] 



> A(n— m)— sup [A(n— m) — A(n— m— x)] (4) 

n— m>0 

Let k = n — m. Eq.© can be written as 

a + (n) — a(m) > \{n — m) — sup[A(/c) — A(/c — x)] 

fc>0 



from which we obtain 

a(n) > a(m) + A(n — m) — y 



(5) 



because e — > and y = sup fc>0 [A(/c) — A(/c — x)]. Since 
Eq.© holds for VO < rre < re, we have 

a(ra) > sup [a(m) + A(n — m) — y] = a<S>A(n) — y. 

0<m<n 



Example 1. Suppose the number of cumulative arrival 
packets of a flow, A(t), is upper-bounded by a(t) + x 
for t > 0, where a(i) = p-t + a. Let n = a(ij. We can 
get the inverse function of a(t), A(n) — ( n ~°") 



We can 



use Eq.© to get y. For V7c > 0, we have 

a) + (k — a — x) 



sup\ — 

fc>0 



< 




P P 

k > a + x, 

a < k < a + x, 

k < a. 



from which we get y = -. Then, we know that for any 



packet, its arrival time satisfies 

(n 



a(n) > sup [a(m) + 

0<m<n 



P 



x 
P 



If o(n) is lower-bounded with respect to some func- 
tion A (re) G Q , we have the following lemma. 

Lemma 3. For function A (re) £ £/, f/iere holds: 

1) the following statements are equivalent: 

a) VO < m < n, a(n) — a{m) > A(n — m) — y 
for \/y > 0; 

b) Vn > 0, a(n) > a(§>A(n) - y for Vy > 0; 

2) if Vn, y > 0, a(n) > a®A(n) — y holds, then we 
have A(t) < A® a(t) + x, where a(t) € Q is the 
inverse function of A (re) and defined as follows 

a(t) = sup{k : X(k) < t} (6) 

and 



x = sup[a(r + y) — a(r) + 1]. 

r>0 



(V) 



Proof: (1) The (a) — > (b) part has been proved in 
Lemma Ht2). We only prove the (b) — ► (a) part. From 
the condition, we have 

a(re) — sup {a(rre) + A(n — m)} + y > 

0<m<n 

which implies 

inf {a(n) — a(rre) — A(n — m)} + y > 0. 

0<m<n 



Thus there holds a(re) — a(rre) > A(n — rre) — y for 
V0 < m < re and Vy > 0. 

(2) For V0 < s < t, we can find rre, re > according 
to the following functions 

A(t) = re = sup{/c : a{k) < t} 

A(s) =m = sup{A; : a(k) < s}. 

Thus, we have A(s, t) = n — m and a(n) — a(m + 1) < 
t — s. From (1), we know that a(n) > a(£)\(n) — y is 
equivalent to a(n) — a(m) > A(n — m) — y. Then we 
have 

t — s > a{n) — a{m + 1) > A(n — rre — 1) — y. 
Taking the inverse function of A (re — m — 1) yields 

re — m — 1 < a(i — s + y) 
Because ,4(s, t) = n — m, we have 
,4(s,t) < a(t-s + y) + l 

= a(t — s) + [a(t — s + y) — a(t — s) + 1] 
< a(t — s) + sup [a(t — s + y) — a(i — s) + 1] 

t-s>0 

Let t = t — s. The above inequality is written as 

A(s, t) < a(t — s) + sup[a(r + y) — a(r) + 1] 

r>0 

= a(t — s) + x 

Since ^.(s, - a(t - s) - x < holds for V0 < s < t, 
we have 

sup [A(s, t) -a(t- s) -x] < 

0<s<t 

from which we further obtain 

A(t) - inf [A(s) + a(t - s)] - x < 0. 

0<s<< 

We then conclude A(t) < A (8) a(t) + x. ■ 

D. Related Space-domain Results 

This sub-section reviews some related space-domain 
results under min-plus algebra O. It is worth high- 
lighting that the following results are for discrete time 
systems with unit discretization step. 

The virtual-backlog-centric (v.b.c) stochastic arrival 
curve model lfl4l is denned based on a probabilis- 
tic upper-bound on cumulative arrival. To ease later 
analysis, the definition of v.b.c stochastic arrival curve 
model presented in this paper is based on the number 
of cumulative arrival packets while not the amount of 



cumulative arrival (in bits) which has been widely used 
in the network calculus literature. 

The v.b.c stochastic arrival curve model explores the 
virtual backlog property of deterministic arrival curve, 
which is that the queue length of a virtual single server 
queue (SSQ) fed with the same flow with a deterministic 
arrival curve is upper-bounded. 

For a flow having arrival curve a(t), we construct a 
virtual SSQ system fed with the same flow. The SSQ 
system has infinite buffer space and the buffer is initially 
empty. Suppose the virtual SSQ system provides service 
a(t) to the flow for all t > 0. Then the unfinished work 
or backlog in the virtual SSQ system by time t is B(t) = 
A(t)—A*(t). The Lindely equation can be used to derive 
B(t), which is 

B(t) = max{0,B(t-l)+A(t-l,t)-a(t-t+l)} (8) 

Eq.© means that the amount of traffic backlogged in the 
system by time t equals the amount of traffic backlogged 
by time t — 1 plus the amount of traffic having arrived 
between t — 1 and t minus the amount of traffic having 
been served between t — 1 and t. By applying Eq.® 
iteratively to its right-hand side, it becomes 

B(t) = sup [A(s,t) -a(t-s)]. (9) 

0<s<t 

If the flow is constrained by arrival curve a(t) + x for 
all t > 0, it follows from Eq.© that the system backlog 
is also upper-bounded by x. The v.b.c stochastic arrival 
curve is defined based on the virtual backlog property. 

Definition 1. (v.b.c Stochastic Arrival Curve). 

A flow is said to have a virtual-backlog-centric (v.b.c) 
sto- 
chastic arrival curve a(t) € Q with bounding function 
f(x) £ G, denoted by A (a, /), if for all t > and 
all x > 0, there holds 

P{ sup [A(s,t) - a(t - a)] > x) < f(x). (10) 

0<s<t 

Eq.dTUb can also be written as follows: 

P{A(t) > A®a(t) +x} < f(x). (11) 

Based on the existing space-domain traffic and server 
models, a lot of results have been derived for stochastic 
network calculus which include the five basic properties 
|4ll as introduced in Sec. HI In this paper, the following 
result is specifically made use of in later analysis and 
hence listed: 

Lemma 4. (Superposition Property). Consider N flows 
with arrival processes Ai(t), i=l,...,N, respectively. Let 



A(t) denote the aggregate arrival process. IfVi, A% ^ v \, 
(ai, fi), then A ^ v b (a, f) with a(t) = Ya=i and 
fix) = h ® • • • ® In(x). 

III. Time-Domain Models 

This section reviews the deterministic arrival curve 
and the deterministic service curve models defined in 
the time-domain. We generalize the deterministic models 
and define time-domain stochastic arrival curve and 
stochastic service curve models. 

A. Deterministic Arrival Curve 

Consider a flow of which packets arrive to a system 
at time a(n). In order to deterministically guarantee a 
certain level of quality of service (QoS) to this flow, 
the traffic sent by this flow must be constrained. The 
deterministic network calculus traffic model in the time- 
domain characterizes packet inter-arrival time using a 
lower-bound function, called arrival curve in this paper 
and defined as follows lfT31 : 

Definition 2. (Arrival Curve). A flow is said to have 
a (deterministic) arrival curve A(n) € Q, if its arrival 
process a(n) satisfies, for all < m < n, 

a(n) — a(m) > A(n — m). (12) 

The arrival curve model has the following triplicity 
principle which will be used as the basis in defining 
the stochastic arrival curve models in the subsequent 
subsections. 

Lemma 5. The following statements are equivalent: 

1) V0 < m < n, a(n) — a(m) > [X(n — m) — x] + 
for Vx > 0; 

2) Vn > 0, sup 0<m<n | [A(n — m) — x\ + — [a(n) — 
a(m)]\ < Ofor \fx > 0; 

3) Vn > 0, sup < m <„ sup < g < m | [A(m - q) - x] + - 
[a(m) - a(q)]} < Ofor Vx > 0, 

where A £ Q. 

Proof: It is trivially true that 

X(n—m) — [a(n)—a(m)} < sup {A(n— m) — [a(n)— a(m)]} 

0<m<n 

from which, (2) implies (1). In addition 

sup {X(n — m) — [a(n) — a(m)]} 

0<m<n 

< sup sup { X(k — m) — [a(k) — a(m)]} 

0<m<n m<k<n 



sup sup {X(k - 

0<KnO<m<fc 

sup sup {A(m 

0<m<n 0<q<m 



m) — [a(k) — a(m)]} 
-q) - [a{m) - a(q)]} 



with which, (3) implies (2). 

For(l)— >(2), it holds since a(n) — a(m) > \{n — m) 
x for VO < m < n. For (2)— >(3), 

sup sup {X(m—q) — [a(m)—a(q)]} < sup [x] = 

0<m<n 0<q<m 0<m<n 



X. 



Thus (1), (2) and (3) are equivalent. ■ 
From Definition |2j the right-hand side of a(n) — 
a(m) > \{n — m)—xm Lemma[5](l) defines an arrival 
curve X(n—m)—x. In addition, we can construct a virtual 
single server queue (SSQ) system that is initially empty, 
fed with the same traffic flow, and has a service curve A 
which makes d(n) < a(g>A(n) (see Definition [6]). Then, 
the delay in the virtual SSQ system is upper-bounded by 
d(n)—a(n) < sup < m<n [A(n— m)— (a(n)— a(m))] < x, 
and the maximum system delay for the first n packets is 
upper-bounded by 



sup {d(m) 

0<m<n 



a(m)} 



< sup sup {A(m 

0<m<n 0<q<m 



q) - [a(m) - a(q)]} < x. 



Example 2. The Generic Cell Rate Algorithm 
(GCRA) 021 with parameter (T, r) is a parallel algo- 
rithm to the Leaky Bucket algorithm and has been used 
in fixed-length packet networks such as Asynchronous 
Transfer Mode (ATM) networks. The GCRA measures 
cell rate at a specified time scale and assumes that cells 
will have a minimum interval between them. Here, T 
denotes the assumed minimum interval between cells 
and r denotes the maximum acceptable excursion that 
quantifies how early cells may arrive with respect to T. It 
can be verified that if a flow is GCRA(T, T)-constrained, 
it has an arrival curve 

A(n) = (T • n — t) + . 

B. Inter-arrival-time Stochastic Arrival Curve 

Lemma |5](1) defines a deterministic arrival curve 
A(n) — x which lower-bounds the inter-arrival time 
between any two packets. Based on this, we define its 
probabilistic counterpart as follows: 

Definition 3. (La.t Stochastic Arrival Curve). A flow 
is said to have an inter-arrival-time (i.a.t) stochastic 
arrival curve A € Q with bounding function h £ Q, 



denoted by a(n) ~# (A, h), if for all < m < n and 
all x > 0, there holds 

pjA(n — m) — [a(n) — a(m)} > x\ < h(x). (13) 

Example 3. Consider a flow with fixed unit packet 
size. Suppose its packet inter-arrival times follow an 
exponential distribution with mean - p . Then, the packet 
arrival time has an Erlang distribution with parameter 
(n,p) ifTTll . And, for any two packets p m and p n , their 
inter-arrival time a(n) — a(m) satisfies, for Vx > 0, 



pj— (n — m) — [a(n) — a(m)] > cc| 



< 1 



n— m— 1 

E 

k=0 



-pv 



{py) k 



-py 



{py) 



n—m—1 



(n 



m 



IV. 



where y = Mn — m) — x. 

The i.a.t stochastic arrival curve is intuitively simple, 
but it has limited use if no additional constraint is 
enforced. Let us consider a simple example to understand 
this problem. Consider a single node with constant per 
packet service time T and its input flow F satisfying 
a(n) ~it (r • n, h) where r > T. Suppose we are inter- 
ested in the delay D(n), where, by definition, D(n) = 
d(n) — a{n). Because the node has constant per packet 
service time T, it has a (deterministic) service curve T ■ n 
which implies d(n) = sup 0<m<r Ja(m) + T ■ (n — m)\. 
Then we have 

D(n) = sup {a(m) + T ■ (n — m)} — a{n) 

0<m<n 

= sup (o(m) + T ■ (n — m) — a(n)} 

0<m<n 



< sup {r • (n 

0<m<n 



m) — [a(n) — a(m)]} (14) 



From Eq. (fT4l) . we have difficulty in further deriving more 
results if no additional constraint is added because we 
only know P{r • (n — m) — [a(n) — a(m)] > x} < 
h(x). When investigating the performance metrics such 
as delay bound and backlog bound in Section IIV-AI we 
meet the similar difficulty. 

C. Virtual-system-delay Stochastic Arrival Curve 

The previous subsection stated the difficulty of ap- 
plying i.a.t stochastic arrival curve to service guarantee 
analysis. This subsection introduces another stochastic 
arrival curve model that can help avoid such difficulty. 
This model is called virtual-system-delay (v.s.d) stochas- 
tic arrival curve. The model explores the virtual system 
delay property of deterministic arrival curve as implied 



by Lemma 0(2), which is that the amount of time a 
packet spends in a virtual SSQ fed with the same flow 
with a deterministic arrival curve is lower-bounded. 

For a flow having deterministic arrival curve, we 
construct a virtual SSQ system fed with the flow, which 
has infinite buffer space and the buffer is initially empty. 
Suppose the virtual SSQ system provides a deterministic 
service curve A to the flow or d(n) = a<g>\(n) for 
all n > 0. The amount of time packet n spends in 
the virtual SSQ system is W s (n) = d(n) — a(n) = 
sup < m < n {A(n — m) — [a(n) — o(m)]}. If the flow is 
constrained by arrival curve X(n) — x for all n > 0, W s 
is also lower-bounded by x. 

Based on the virtual system time property, we define 
virtual-system-delay (v.s.d) stochastic arrival curve to 
characterize the arrival process. 

Definition 4. (v.s.d Stochastic Arrival Curve). A flow 
is said to have a virtual-system-delay (v.s.d) stochastic 
arrival curve A G Q with hounding function h G Q, 
denoted by a(n) (A, h), if for all < m < n and 
all x > 0, there holds 



P< sup {A(n — m) — [a(n) — a(m)\] > x > < h(x). 

0<m<n ' 

(15) 



(16) 



Eq.([15l) can also be written as 



P{a®\(n) - a(n) > x) < h(x). 



a<g>A(n) can be considered as the expected time that 
the packet would arrive to the system if the flow had 
passed through the virtual SSQ with service curve A(n). 
x denotes the difference between the expected arrival 
time and the actual arrival time. Eq.(fT6b characterizes 
this difference x by introducing a bounding function 
h(x). 

Example 4. Consider a flow with the same fixed 
packet size. Suppose all packet inter-arrival times are 
exponentially distributed with mean ^. Based on the 
steady-state probability mass function (PMF) of the 
queue-waiting time for an M/D/l queue [18], we say 
that the flow has a v.s.d stochastic arrival curve a(n) ~ vd 
(D ■ n, h exp ) for V£> < i , with p = fi- D and 

h e *P(x) = 1 - (1 - p) ]T e-rt- x > m jJ 

i=Q l - 

where, [x/D\ denotes the greatest integer less than or 
equal to x/D. 



The following theorem establishes relationships be- 
tween i.a.t stochastic arrival curve and v.s.d stochastic 
arrival curve. 

Theorem 1. 1) If a flow has a v.s.d stochastic arrival 
curve A € Q with hounding function h G Q, then 
the flow has an i.a.t stochastic arrival curve A € Q 
with the same bounding function h € Q. 
2) Inversely, if a flow has an i.a.t stochastic arrival 
curve A G Q with bounding function h € T, it also 
has a v.s.d stochastic arrival curve A_^ G Q with 
bounding function h r ' G Q where 



A_ r/ (n) = A(n) — r/ ■ n 
hJ'ix) = \h(x) + - f°° h(y)dy 
for Vt? > 0. 

Proof: The first part follows from that 

X(n—m) — [a(n)—a(m)} < sup {A(n— m) — [a(n)— a(m)]} 

0<m<n 

holds for V0 < m < n. For the second part, there holds 

sup {A_ r? (n — m) — [a(n) — a(m)]} 

0<m<n 

< st sup {\- n {n - m) - [a(n) - a(m)]} + 

0<m<n 

Since for Vx > 0, 

P{{X(n — m) — r\ ■ (n — m) — [a(n) — a(m)]} + > x} 
= P{{X(n — m) — rj-(n — m) — [a(n) — a(m)]} > x} 
< h(x + rj ■ (n — m)), 

we have 

P\ sup {\- v (n — m) — [a(n) — a(m)]} > x > 

0<m<n ' 



< 



p|{A_^(n — m) — [a(n) — a(m)]} + > x| 

m=0 

n n 

< \_. h{x + rj ■ (n — m)) = 2_, + V ' ^) 



m=0 

oo 



k=0 

oo 



< h(x + r] ■ k) = h(x) + h(x + r\ ■ k) 



k=0 



k=l 



l r 00 

< h(x) + - / h{y)dy. 

V Jx 



(17) 



which is meaningful only when Eq.(fT71) is upper- 
bounded by one. The 1-fold integration of h(x) is 
bounded by one because the condition assumes h G T as 



for the HI. Then the second part follows from Eq.dTTT). 

■ 

Note that in the second part of the above theorem, 
h(x) G T while not G Q. If the requirement on the 
bounding function is relaxed to h(x) G G, the above 
relationship may not hold in general. 

The v.s.d stochastic arrival curve has a counterpart 
defined in the space-domain, the v.b.c stochastic arrival 
curve as defined in Definition Q] The following theorem 
establishes relationships between these two models. 

Theorem 2. 1) If a flow has a v.b.c stochastic arrival 
curve at(t) G Q with bounding function f(x) G Q, 
the flow has a v.s.d stochastic arrival curve A(n) G 
Q with bounding function h(y) G Q, where A(n) = 
inf{r : a(r) > n} and h(y) = /(sup r > [a(r + 
y) - a(r) + 1]). 
2) If a flow has a v.s.d stochastic arrival curve 
A(n) G Q with bounding function h(y) G Q, the 
flow has a v.b.c stochastic arrival curve a(t) G Q 
with bounding function f{x) G Q, where a(t) = 
sup{fc : \{k) < i} and f(x) = h(sup k>0 [X(k) — 
A(fc -*)])■ 

Proof: (1) From Lemma |2l we know that for Vx, t > 
0, event {A(t) < A (g> a(t) + x} implies event {a(n) > 
a®A(n) — y} where y is obtained from Eq.©. Thus, 
there holds 

P{A{t) < A®a(t)+x} < P{a{n) > a®\(n)-y}. 

We further have 

P{A(t) > A®a(t) + x} > P{a(n) < a®\(n)-y}. 

From the condition that the flow has a v.b.c stochastic 
arrival curve a(t), we know P{A(t) > A®a>(t) + x} < 
f(x). According to Eq.©, we obtain 

P{a(n) < a®\(n)-y} < /(sup[a(r+y)-a;(r) + l]). 

T>0 

(2) From Lemma [51 we know that for Vy > 0, event 
{a(n) > a(&\(n) — y} implies event {A(t) < A®a(t) + 
x} where x is obtained from Eq.©. Thus, there holds 

P{a{n) > a®\{n)-y} < P{A{t) < A®a(t) + x} 

We further have 

P{a(n) < a®\(n)-y} > P{A{t) > A<g>a(t)+x} 

From the condition that the flow has a v.s.d stochastic 
arrival curve A(n), we know P{a(n) < a(§\(n) — y} < 
h(y). According to Eq.((3j, we have 

P{A(t) > A®a(t) + x} < h(sup[X(k) - X(k - x)]) 

k>0 

and complete the proof. ■ 



D. Maximum-(virtual)-sy stem-delay Stochastic Arrival 
Curve 

The maximum-(virtual)-system-delay (m.s.d) stochas- 
tic arrival curve explores the maximum virtual system 
delay property of deterministic arrival curve implied 
by Lemma |5](3), which is that the maximum system 
delay of a virtual SSQ fed with the same flow with a 
deterministic arrival curve is lower-bounded. 

Similar to the discussion for v.s.d stochastic arrival 
curve, for a flow having arrival curve, we construct a 
virtual SSQ system fed with the flow, which has infinite 
buffer space and the buffer is initially empty. Suppose 
the virtual SSQ system provides a deterministic service 
curve A to the flow or d(n) = a<§)A(n) for all n > 0. 
The maximum system delay in the virtual SSQ system 
for the first n arrival packets as sup 0<m<n W s (m) = 
sup < m <n sup < g < m {A(m - q) - [a{mj- a(q)}}. If the 
flow is constrained by arrival curve A(n) — x for all 
n > 0, the maximum system delay in the virtual SSQ is 
also upper-bounded by x. 

Based on the maximum virtual system delay property, 
we define m.s.d stochastic arrival curve model. 

Definition 5. (m.s.d Stochastic Arrival Curve). 

A flow is said to have a maximum-(virtual)- system- 
delay (m.s.d) stochastic arrival curve A(n) G Q with 
bounding function h(x) G Q, denoted by a(n) ^ m d 
(A, h), if for all < m < n and all x > 0, there holds 

P< sup sup {X(m—q) — [a(m)—a(q)]}>x><h(x). 

0<m<n 0<q<m 

(18) 

E. Deterministic Service Curve 

To provide service guarantees to an arrival-constrained 
flow F, the system usually needs to allocate a minimum 
service rate to F. A guaranteed minimum service rate 
is equivalent to a guaranteed maximum service time for 
each packet of the flow, and accordingly the packet's 
departure time from the system is bounded. Because 
packets of the same flow are served in FIFO manner, 
any packet p n from this flow will depart by d(n) which 
is iteratively defined by 

d(n) = max[a(n), d(n — 1)] + 5(n) (19) 

with <i(0) = 0, where 5(n) is the service time guaranteed 
to p n . By applying Eq.(fT9l iteratively to its right-hand 
side, it becomes 

n 

d(n) = sup [a(m) + ^ S(i)} (20) 

0<m<n • 



where J2i= m ^(*) * s tne guaranteed cumulative service 
time for packet p m to p n . Suppose we can use a 
function "y(n — m) to denote J2?= m <K*)> ie - 7( n — m ) = 
X)r=m^W- Then, Equation(f20| becomes 

d(n) = sup [a(m) + 7(72 — m)] = a®^y(n) 

0<m<n 

which provides a basis for the following time-domain 
(deterministic) server model that charaterizes the service 
using an upper bound on the cumulative packet service 
time ifBl : 

Definition 6. (Service Curve). Consider a system S with 
input process a(n) and output process d{n). The system 
is said to provide to the input a (deterministic) service 
curve j(n) E Q, if for Vn > 0, 

d(n) < a®j(n). (21) 

The (deterministic) service curve model has the fol- 
lowing duality principle: 

Lemma 6. For Vx > 0, d(n) — a<S>7(n) < x for all 
n > 0, if and only if sup 0<m<n [d(n) — 0(8)7(71)] < x for 
Vn > 0, where 7 £ 

Proof: For the "if part, it holds because d(n) — 
a®7(n) < sup 0<m<Tl [d(n) — a®j(n)]. For the "only if 
part, from d(n) — a®^(n) < x for Vn > 0, we have 
sup < m <n[d(n) - o®7(n)] < sup < m < n [x] = x. ■ 
By the definition of service curve, the first part of 
Lemma [6] defines a service curve j(n) + x. Lemma [6] 
states that if a server provides service curve 7(n) + x, 
then sup 0<m<n [d(m) — a®7(m)] < x holds, and vice 
versa. In this sense, we call Lemma [6] the duality 
principle of service curve. 

F. Stochastic Service Curve 

For networks providing stochastic service guarantees, 
following the principle of Eq.(l20l). we have the following 
expression for the expected departure time of packet p n 

n 

d(n) = sup [a(m) + ^ (5(i) + e(t))] 
o<m<n i=m 

where we assume 5(i) is the deterministic part while e(i) 
the random part in the total service time 5(i)+e(i) guar- 
anteed to packet p l . We call e(i) stochastic error term 
associated to 6(i). Here, e(n) is introduced to represent 
the additional delay of p n due to some randomness. For 
example, an error-prone wireless link is often considered 
to operate in two states. If the link is in 'good' condition, 
it can send and receive data correctly; if the link is in 
'bad' condition due to errors, the data that should be sent 



immediately has to be queued longer until the channel 
changes to 'good' condition. Then, e(n) in this case 
represents the time period in which the channel is in 
'bad' condition between the time when p n_1 has been 
sent correctly and the time when p n can be sent. 

With the consideration of the stochastic error term, 
the (deterministic) service curve can be extended to a 
stochastic version as follows: 

Definition 7. (i.d Stochastic Service Curve). 

A system is said to provide an inter-departure time (i.d) 
stochastic service curve 7 G Q with hounding function 
j € G, denoted by S if for all n > and all 

x > 0, there holds 

P\d(n) - a®7(n) > x\ < j(x). (22) 

Example 5. Consider two nodes, the sender and the 
receiver, communicate through an error-prone wireless 
link. Packets have fixed-length. Packets arriving to the 
sender node are served in FIFO manner. Assume the 
guaranteed per-packet service time is 5 without any 
error. To simplify the analysis, assume the time slot 
length equals S. The sender sends packets correctly only 
when the link is in 'good' condition. If the link is 
in 'bad' condition, no packets can be sent correctly. In 
addition, the sender can send the head-of-queue packet 
only at the beginning of a time slot, i.e., the time period 
during which the link is in 'bad'condition should be 
an integer times of 5. The probability that a packet 
can be sent correctly is determined by packet error rate 
(PER). PER is determined by the packet length and the 
bit error rate (BER). Here, we assume packet errors 
happen independently and the same PER denoted by 
P e is applied to all packets. The successful transmission 
probability of one packet is hence 1 — P e . 

Suppose P{A(n) = *} = P e i_1 (l - -Pe), i > 1, where 
A(n) represents the number of time slots necessary to 
successfully send the n th packet with respect to the 
successful transmission probability 1 — P e . The number 
of time slots necessary to successfully send n packets is 
J2k=i A(fc) which has the negative binomial distribution 

p{± m=i }J (i: 1 i)( i - , w"^ 

k=1 { 0, i < 

Then the sender provides to its input a stochastic service 
curve 7 which has the following distribution 



where r is the guaranteed service time to successfully 
send n packets and \x] denotes the smallest integer 
greater than or equal to x. 

We can find no < n such that a<g>7(n) takes its 
maximum value, i.e., a®j(n) = a(no) + j(n — no + 1). 
From Eq. (|22l) . we have 



P{7(n — no + 1) < d(n 
where 

r d(n) — a(Ti ) — x -| ^ 



a(n ) - < 



n — 1 



(1 - P e )"P ( 



In Sec. [TV] we show that many results can be derived 
from the i.d stochastic service curve model. However, 
without additional constraints, we have difficulty in prov- 
ing the concatenation property for i.d stochastic service 
curve. To address this difficulty, we introduce a stronger 
definition in the following subsection. 

G. Constrained Stochastic Service Curve 

The constrained stochastic service curve model is 
generalized from the (deterministic) service curve model 
based on its duality principle. From Lemma [6l we know 
that a system with input a(n) and output d(n) has a 
service curve 7(n) if and only if for Vn > 0, 



sup {d(m) — a®7(m)} < x. 

0<m<n 



(23) 



Inequality (1231 provides the basis to generalize the 
(deterministic) service curve model to the constrained 
stochastic service curve defined as follows: 

Definition 8. (Constrained Stochastic Service Curve). A 

system is said to provide a constrained stochastic service 
curve (c.s) 7 G Q with bounding function j G Q, denoted 
by S ~ cs (7, j), if for Vn, x > 0, there holds 

P< sup [d(m) — a<8>7(m)] > x > < j(x). (24) 

0<m<n J 

The following theorem establishes relationships be- 
tween i.d stochastic service curve and c.s stochastic 
service curve. 

Theorem 3. I) If a server S provides to its input 
a{n) a c.s stochastic service curve 7(n) with 
bounding function j{x) G Q, it provides to the 
input a(n) an i.d stochastic service curve 7(n) 
with the same bounding function j(x) G Q, i.e., 

S ~id (1,3); 

2) If a server S provides to its input a(n) an i.d 
stochastic service curve 7(n) with bounding func- 
tion j{x) G T, it provides to the input a{n) a c.s 



j{y)dy 



stochastic service curve 7 +r; (n) = ^(nj+rj-n with 
bounding function ffe) G T for V77 > 0, where 

-y J x—rj-n 

Proof: The first part follows since there always 
holds d(n) — a®^(n) < sup < m < n {d(m) — 0(87(771)}. 
For the second part, there holds for V0 < m < n, 

a(X>7 +r? (m) > a(g>7(m) + y ■ m — y ■ n 

and then 

d(m) — a®j +v (m) < d{m) —a®'y(rn) — y-m + y-n. 
Thus, we obtain 

P{ sup {d(m) — a<g)7 +J) (m)} > x} 

0<m<n 

< P{ sup \d{m) — a®~f{m) — y(m)] + > x — y ■ n} 

l<m<n 

for which when x — y ■ n < 0, the right hand side is 
equal to 1. In the following, we assume x — y ■ n < 
under which, there holds 

P{ sup {d{m) — a(g>7 +?? (m)} > x) 

0<m<n 
n 

< E P{ [d(m) — a(§>7(m) — y(m)] > x — y ■ n} 



m=l 

n 



" 1 f n 

< E j( x — V ■ n + y ■ m) < — / 

m=l ^ Jx- V -n 



j{y)dy. 



Because the probability is always not greater than 1, the 
second part follows from the above inequality. ■ 
In the second part of the above theorem, j(x) G T 
while not G Q. If the requirement on the bounding 
function is relaxed to j(x) G Q, the above relationship 
may not hold in general. 

IV. Basic Properties 

This section presents results derived from the time- 
domain traffic models and server models introduced 
in Sec. [TlTJ Particularly, we investigate the five basic 
properties introduced in Sec. H which are service guar- 
antees including delay bound and backlog bound, output 
characterization, concatenation property and superposi- 
tion property. However, some properties can directly be 
proved only for the combination of a specific traffic 
model and a specific server model. This explains why 
we need to establish the various relationships between 
models in Sec. [Till With these relationships, we can 
obtain the corresponding results for models which we 
are interested in. 



A. Service Guarantees 

This subsection investigates probabilistic bounds on 
delay and backlog under the combination of v.s.d 
stochastic arrival curve and i.d stochastic service curve. 

We start with deriving the bound on delay that a packet 
would experience in a system. 

Theorem 4. (Delay Bound). Consider a system S 
providing an i.d stochastic service curve 7 G Q with 
bounding function j G Q to the input which has a v.s.d 
arrival curve A G Q with hounding function h G Q. Let 
D{n) = d{n) — a{n) be the delay in the system of the 
n th (> 0) packet. For Vx > 0, D(n) is bounded by 

P{D{n) > x} <j®h(x -70 A(0)). (25) 

Proof: For Vn > 0, there holds 

d(n)—a(n) = [d(n)—a®j(n)] + [a®j(n)—a(n)] 
= \d[n) — a®7(n)] + sup {X(n—m) — [a(n) — a(m)] 

0<m<n 

+ 7(71. — m) — A(n — m)} 

< [d(n)— a®7(n)] + sup {A(n— m)— [a(n)— a(m)] } 

0<m<n 

+ sup {7(n — m) — A(n — m)} 

0<m<n 

< [d(n) — a<8>7(n)]+ sup {A(n — m)— [a(n) — a(m)] } 

0<m<n 

+ sup (7(/c) - X(k)}. (26) 

fc>0 

The right-hand side of Eq.(l26b implies a sufficient 
condition to obtain P{D(n) > x}, which is that 

P{d(n) — a(§>7(n) > x} and p| sup 0<m<n {A(n — m) — 

[a(n) — a(m)] } > x| are known. To ensure the system's 
stability, we should also have 

lim ~[j(k) - \(k)] < 0. (27) 

In the rest of the paper, without explicitly stating, we 
shall assume inequality (|27T ) holds. From Lemma Q] and 

sup fc > {l(k) 

— X(k)} =70 A(0), we conclude 

P{D(n) > x} < j h(x - 7 A(0)). 

■ 

Next, we consider backlog bound of a system. By 
definition, the backlog in the system at time t > is 
B(t) = A[t) - A*{t). If a{n) is the arrival time of the 



latest packet arriving to the system by time t, then B(t) 
is 

B(t) < inf {A; > : d(n — k) < a(n)}. (28) 

Eq.(f28l) implies that, for Vx > 0, if B(t) > x, there 
must be a(n) < d{n—x). Thus event {B{t) > x} implies 
event {a(n) < d(n-x)} and P{B(t) > x} < P{a(n) < 
d(n—x)}. Then we have the following result for backlog. 

Theorem 5. (Backlog Bound). Consider a system S 
providing an i.d stochastic service curve 7 G Q with 
bounding function j G Q to the input which has a v.s.d 
stochastic arrival curve A G Q with bounding function 
h G (y. The backlog at time t (t > 0), B(t), is bounded 
by 

P{B(t)> tf(A,7 + x)} <j®h{x) (29) 

for any x > 0, where, H(X,j + x) = sup n>0 { inf[k > 
: 7(72 — k) + x < A(n)]} is the maximum horizontal 
distance between functions A(n) and j(n) +x for Vx > 
0. 

Proof: Similar to prove the delay bound, we have 

d(n—x)—a(n) = [d(n— x)— a07(n— x)]+[a07(n— x)— a(n)] 
= [d{n— x)— a(§>7(n— x)] + sup {X(n—k) — [a(n)—a(k)] 

0<k<n—x 

+ 7(72 — x — k) — X(n — k)} 
< \d(n— x)— a(§>7(n— x)]+ sup {X(n—k)—[a(n+x) 

0<k<n—x 

—a(k)]}+ sup {j(n — x — k) — X(n — k)} 

0<k<n—x 

Let v = n — k. The above inequality is written as 

d(n—x)—a(n) < [d(n— x)— a07(n— x)] + sup {X(n—k) 

0<k<n 

— [a(n) — a(k)] }+ sup {j(v — x) — X(v)} 

x<v<n 

Let x = H(X, 7 + y), we have 

d(n-h(X,^+y))-a(n) < [d(n—h(X,j+y))- 

a07(n-/i(A,7+y))]+ sup {X(n-k)-[a(n)-a(k)]}-y 

0<k<n 

(30) 

Under the same conditions as analyzing the delay, we 
obtain 

P{B(t) > H(X, 7 + x)} < j h{x). 



B. Output Characterization 

This subsection presents the result for characterizing 
the departure process from a system. 

Theorem 6. (Output Characterization). Consider a sys- 
tem S provides an id stochastic service curve 7(71) G Q 
with bounding function j(x) G Q to its input which 
has a v.s.d stochastic arrival curve A(n) G Q with 
bounding function h(x) G Q. The output has an i.a.t 
stochastic arrival curve X<2)j(n — m) with bounding 
function j <g> h{x) G Q. 

Proof: For any two departure packets m < n, there 
holds 

d(n)—d(m) > a{n)—a®^f{m)+a® , y(m)—d{m) 

— [d(n) — d(m)] < [d(m) — a<g>7(m)] +a<g>7(m)— a(n) 

< [d(m)— a®7(m)] + sup {X(n—k) — [a(n)—a(k)]} 

0<k<n 

+ sup {7(f) — A(n — m + v)} 

0<v<m 

= \d(m)— a®7(m)] + sup {X(n—k) — [a(n)—a(k)]} 

0<k<n 

— inf {A(n — m + v) — 7(f)} 

0<v<m 

Adding info<„< m {A(n — m + v) — j(v)} to both sides 
of the above inequality, we get 

inf {\(n — m+v)—j(v)} — [d(n) — d(m)] 

II' v 

< [d(m) — o®7(m)] + sup {A(n — k) — [a(n) — a(k)\] 

0<k<n 

With the same conditions as analyzing delay, we con- 
clude 

-PjA07(n — m) — [d(n) — d(m)\ > x| < j (gi h(x) 



C. Concatenation Property 

The concatenation property uses an equivalent system 
to represent a system of multiple servers connected in 
tandem, each of which provides stochastic service curve 
to the input. Then the equivalent system provides the 
input a stochastic service curve, which is derived from 
the stochastic service curve provided by all involved 
individual servers. 

Theorem 7. (Concatenation Property). Consider a flow 
passing through a network of N systems in tandem. If 
each system k(= 1,2, ...,N) provides a c.s stochastic 
service curve S k ~ cs ('J k ,j k ) to its input, then the 



network guarantees to the flow a c.s stochastic service 
curve S ~ cs (7, j) with 

7(71) = 7 1 <g>7 2 (g> • ■ • ®7 N (n) 



j(x) = j 1 (g> j 2 



(31) 



(32) 



Proof: We shall only prove the two-node case, from 
which, the proof can be easily extended to the iV-node 
case. The departure of the first node is the arrival to 
the second node, so d}{n) = a 2 {n). In addition, the 
arrival to the network is the arrival to the first node, i.e., 
a(n) = a 1 (n), and the departure from the network is 
the departure from the second node, i.e., d(n) = d 2 (n), 
where, a(n) and d(n) denote the arrival process to and 
departure process from the network, respectively. We 
then have, 



sup {d(m) — a(8)7 1 ®7 2 (m)} 

0<m<n 

= sup {d 2 (m) - (a 1 (g>7 1 )(g)7 2 (m)} 

0<m<n 



(33) 



Now let us consider any m, (0 < m < n), for which we 
get, 

d 2 (m) - (a 1 (g)7 1 )«)7 2 (m) 
= d 2 (m)- sup |a 1 (g)7 1 (/c)+7 2 (m-A;)-d 1 (/i;)+a 2 (A:)} 

0<fc<m 

= d 2 {m)+ inf {d 1 (k)-a 1 ^ 1 (k)-j 2 {m-k)-a 2 (k)} 

0<k<m 

< sup {d 1 {k)-a l ^ l {k)}+d 2 {m) 

0<k<m 

+ inf {-[a 2 (k)+ 1 2 {m-k)}} 

0<k<m 



< sup {d (k) — a (g>7 (&)}+ [d (m) — a ®7 (m)] 

0<AKm 

(34) 

Applying Eq.(l33l to Eq.(T34l. we obtain 

sup {d 2 (m)-(a 1 (g)7 1 )(§)7 2 (m)} 

0<m<n 



< sup {d 1 (k)-a l ®j l (k)}+ sup {d^m)-^®^^)} 

0<k<n 0<m<n 

(35) 

with which, since both nodes provide c.s stochastic 
service curve to their input, the theorem follows from 
Lemma Q] and the definition of c.s stochastic service 
curve. ■ 



D. Superposition Property 

The superposition property means that the superposi- 
tion of flows can be represented using the same traffic 
model. With this property, the aggregate of multiple 
individual flows may be viewed as a single aggregate 
flow. Then the service guarantees for the aggregate flow 
can be derived in the same way as for a single flow. 

First, we only consider the aggregate of two flows, F\ 
and F2. Let a\(n), a>2 (n) and o(n) be the arrival process 
of Fi, F2 and the aggregate flow Fa, respectively. 

For any packet p n of the aggregate flow Fa, it is either 
the m th packet from flow F\ or the [n — m) th packet 
from flow F2, where m € [0, n), i.e. 

a(n) = max{ai(m), (i2{n — m)} 

For example, a(l) is either max[ai(0), 02(1)] or 
max[oi(l), 02(0)] and the minimum of these two possi- 
bilities, 

a(l) = inf { max[ai(0), 02(1)], max[ai(l), 02(0)]}. 

We can see another example 

a(2) = inf { max[ai(0), 02(2)], max[ai(l), 02(1)], 

max[ai(2), a 2 (0)]}. 

Essentially, we have for any packet n of the aggregate 
flow 

a(n) = inf { max[ai(m), a2(n — m)]}. (36) 

0<m<n 

We generalize the result to the superposition of iV(> 2) 
flows 

a(n) = inf { max[ai (mi), 02(7^2), 

y~] rrii=n 

N-l 

a N (n - rrii)]}. (37) 
i=i 

From Eq.(l37l). it is difficult to directly characterize the 
packet inter-arrival time of the aggregate flow. We know 
that if a flow has a v.s.d stochastic arrival curve, with 
Theorem 02), this flow has a v.b.c stochastic arrival 
curve, for which the superposition property holds 
Thus, we can indirectly prove that the superposition 
property holds for the v.s.d stochastic arrival curve. 

If flow i has a v.s.d stochastic arrival curve aj(n) ^ vd 
(Xi,hi) i = 1,2, ...,N, from Theorem [32), flow i 
has a v.b.c stochastic arrival curve oti(t) with bounding 
function /j(x) = hi( sup T > [a^r + y) - a^r) + 1]) 
, where aj(i) = sup{A: : Xi(k) < t}. According to 
Lemma |4l the aggregate flow has a v.b.c stochastic 



arrival curve a(t) = Ya=i a i(t) with bounding function 
f(x) = fi gj • • • <g> /at(x). We apply Theorem 01) and 
obtain the following result: 

Theorem 8. Consider N flows with arrival processes 

ai{n) ~ vd (Aj, hi), i = 1, N. For the aggregate of 
these flows, there holds a(n) r^ vd (A, h) with A(n) = 
inf{r : a{r) > n} and 

h(y) = f(snp{X(k)-X(k-x))), 

K k>0 ' 

where a(t) = J2iLi a iif) an d f(%) = /1 ® • • • ® In(x) 
with 

cti(t) = sup{k : Xi(k) < t} 

fi(x) = hi(sup[ai(T + y)- a^r) + 1]). 

T>0 

E. Leftover Service Characterization 

This subsection explores the leftover service char- 
acterization under aggregate scheduling. To ease the 
discussion, we consider the simplest case when there are 
two flows competing resource in a system under FIFO 
aggregation. Suppose that if packets arrive to the system 
simultaneously, they are inserted into the FIFO queue 
randomly. Consider a system fed with a flow Fa which 
is the aggregation of two constituent flows F\ and F^. 
Suppose both the service characterization from the server 
and traffic characterization from F2 are known. We are 
interested in characterizing the service time received by 
F\ , with which per- flow bounds for F% can be then easily 
obtained using earlier results derived in the previous 
subsections. 

Theorem 9. Consider a system S with input Fa that 
is the aggregation of two constituent flows F\ and F2. 
Suppose F2 has a (deterministic) arrival curve X2(n) 6 
Q, and the system provides to the input an i.d stochastic 
service curve 7 € Q with bounding function j(x) G Q. 
Then if 7(n + sup[g : X2(q) < di(n)]) € Q, F\ receives 
an i.d stochastic service curve 7(71 + sup[g : X^iff) < 
ai(n)]) with the same bounding function j(x). 

Proof: Suppose packet p\ is the (n + m) th packet 
of Fa, i.e., a(n + m) = a\(n), where m represents the 
number of packets from F2. As the system provides an 
i.d stochastic service curve 7(71) to the aggregate flow 
Fa, there holds 

P{d(n + to) — a®7(n + to) > x} < j(x). 

oi(n) = a(n + m) indicates 02 (to) < a\(n). Let to = 
sup[(/ : X2{q) < a\[n)]. As A2 is the (deterministic) 



arrival curve of F2, we have rh > m because of a% (m) > 
\2(m). Then 7(71 + m) > y(n + m). Let 71 (n) = 7(71 + 
m). From 71 (n) > 7(72 + m), we have ai<8>7i(n) > 
a<8>7(n + m). As d(n + m) = di{n), there holds 

di(n) — ai<8>7i(n) < d(n + m) — a®y{n + m). 
Thus, we conclude 

P{di(n) - ai<8>7i(n) > x} < 
and complete the proof. ■ 

F. Discussion 

In this section, we have presented the five basic prop- 
erties of stochastic network calculus under various traffic 
models and server models defined in the time-domain 
and introduced some simple applications. For example, a 
GCRA-constrained flow has a deterministic arrival curve. 
If a flow's packet inter-arrival times are exponentially 
distributed, then this flow has a v.s.d stochastic arrival 
curve. The service process of an error-prone wireless link 
can be modeled by an i.d stochastic service curve. 

For each basic property, we investigated one com- 
bination of a specific traffic model and a specific 
server model. Particularly, we proved that the service 
guarantees and the output characterization hold for the 
combination of v.s.d stochastic arrival curve and i.d 
stochastic service curve. For the concatenation property, 
we investigated the case that all servers provide the 
constrained service curve to their input but did not 
specify the type of arrival curve. In order to prove the 
supeiposition property, we used the transformation be- 
tween v.s.d stochastic service curve and v.b.c stochastic 
service curve. The leftover service characterization was 
only proved for the combination of deterministic arrival 
curve and i.d stochastic service curve. 

With the relationships and transformations among 
models established in Sec. [IlTJ these five properties may 
be directly or indirectly proved for other combinations 
of traffic models and server models. For example, it is 
easy to prove the service guarantees and output charac- 
terization for the combination of m.s.d stochastic arrival 
curve and i.d stochastic service curve. Considering space 
limitation, these results are not included. However, to 
prove the concatenation property and superposition prop- 
erty for other server models and traffic models, it will 
require additional transformations among models. For 
the leftover service characterization, we may need more 
constraints or transformations when proving it for other 
combinations of traffic models and server models. We 
leave these as our future work. 



V. Conclusion 

For stochastic service guarantee analysis, we intro- 
duced several time-domain models for traffic and service 
modeling. The essential idea of them is to base the model 
on cumulative packet inter-arrival time for traffic and 
on cumulative service time for service. Simple examples 
have been given to demonstrate the use of them. Based 
on the proposed time-domain models, the five basic 
properties for stochastic network calculus were derived, 
with which, the results can be easily applied to both the 
single-node and the network cases. 

As Example 5 showed, we can directly obtain the 
service curve in the time-domain. With the result of 
service guarantees, the probabilistic delay bound and 
probablistic backlog bound can be readily obtained. We 
believe, the proposed time-domain models and derived 
results can be particularly useful for analyzing stochastic 
service guarantees in systems, where the behavior of 
a server involves some stochastic processes which can 
be directly characterized in the time-domain, while it is 
difficult to characterize such stochastic processes in the 
space-domain. Such systems include wireless links and 
multi-access networks where backoff schemes may be 
employed. 

In this paper, we only analyzed a simple case of 
wireless network to illustrate how to apply the proposed 
server model to characterize the service process of a 
wireless node. The future work is to investigate the 
performance of some typical contention-based multi- 
access networks including IEEE 802.11 networks. 
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